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ABSTRACT 

We present a high-resolution radio survey of the Sloan Digital Sky Survey (SDSS) Southern Equa- 
torial Stripe, a.k.a. Stripe 82. This 1.4 GHz survey was conducted with the Very Large Array (VLA) 
primarily in the A-configuration, with supplemental B-configuration data to increase sensitivity to 
extended structure. The survey has an angular resolution of 1.8" and achieves a median rms noise 
of 52 /iJy beam~^ over 92 deg^. This is the deepest 1.4 GHz survey to achieve this large of an area, 
filling a gap in the phase space between small, deep and large, shallow surveys. It also serves as a 
pilot project for a larger high- resolution survey with the Expanded Very Large Array (EVLA). We 
discuss the technical design of the survey and details of the observations, and we outline our method 
for data reduction. We present a catalog of 17,969 isolated radio components, for an overall source 
density of ^195 sources deg~^. The astrometric accuracy of the data is excellent, with an internal 
check utilizing multiply-observed sources yielding an rms scatter of 0.19" in both right ascension and 
declination. A comparison to the SDSS DR7 Quasar Catalog further confirms that the astrometry 
is well tied to the optical reference frame, with mean offsets of 0.02" ± 0.01" in right ascension, and 
0.01" ± 0.02" in declination. A check of our photometry reveals a small, negative CLEAN-like bias 
on the level of 35 /iJy. We report on the catalog completeness, finding that 97% of FIRST-detected 
quasars are recovered in the new Stripe 82 radio catalog, while faint, extended sources are more likely 
to be resolved out by the resolution bias. We conclude with a discussion of the optical counterparts 
to the catalog sources, including 76 newly-det ected radio quasars. The full catalog aswell as a search 



page and cutout server are available online at http://third.ucllnl.org/cgi-bin/stripe82cutout 
Key words: catalogs - radio continuum: general - surveys 



1. INTRODUCTION 

The history of radio surveys has played out as the 
usual tradeoff between depth and width, with most 
surveys favoring one extreme at the loss of the other. 
On one end, large area radio surveys like Faint Im- 
ages of the Radio Sky at Twen ty centimeters (FIRST; 
iBecker. White, fc Helfanl [l995h and the NRAO VLA 
Skv Survev rNVSS: ICondon et al.l[l998h have been of un- 
questionable importance to the astronomy community. 
These surveys revolutionized the field of radio astron- 
omy with their detailed censuses of the mJy radio pop- 
ulation. Still, they are relatively shallow surveys, with 
typical rms values of 0.15 mJy beam~^ and 0.45 mJy 
bearn"^, respectively. 

At the other end of the spectrum are relatively small 
but deep radio surveys, with pencil beam surveys 
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being the most extreme example. iRichardi (|2OOO0 
surveyed a 1 deg^ region containing the Hubble Deep 
Field North (HDF-N) at 1.4 GHz and achieved a 
completeness limit of 40 ^Jy. Other exa mples of 
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deep 1.4 GHz surveys include those done_ 

(2000); Hopkins ctal. (2003); Condon et al.l 

IBondi et al. (2003); Scvmour ct al (2004); Norris ct aL| 
12005, 2006); Simuson et al. (2006 ): Ivison et al. (2007); 
Schinnerer et al. (2007); Ow en fc MorrisonI (12008): 
Kell erman et al.l (1200811 : iMiller et al.l (|2008D : and 
:Morrison et al.l (|2010D . The "widest", arguably, of 
the deep surv eys, the Phoenix Deep Survey (PDS; 
iHopkins et all [2003). covers 4.5 deg^ and has an rms 
of 12 ^Jy in its most sensitive region. Such deep 
surveys have many important applications, including 
(but not limited to) investigating the dust-unbiased star 
formation in galaxies, the evolution of the radio-FIR 
correlation, and the nature of the somewhat controver- 
sial sub-mjy population; however, they can suffer from 
small number statistics depending on which sources are 
of interest. 
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The future of radio astronomy is bright, or rather, 
faint, as it turns out. ASKAP, the Australian Square 
Kilometer Array Pathfinder, is slated to cover the entire 
Southern sky with 10" resolution to 10 /iJy rms. APER- 
TIF, the Westerbork equivalent of ASKAP, will cover the 
Northern sky with 10" resolution to 20 /xJy rms. The 
Expanded Very Large Array (EVLA), which as already 
begun operations with reduced capabilities, will be 5-20 
times more sensitive than the VLA and give point-source 
sensitivity better than 1 /iJy between 2-40 GHz. The 
Karoo Array Telescope (MeerKat), which South Africa 
is building as a precursor instrument for the SKA, will 
provide a southern hemisphere complement to the EVLA 
from the L-band to the X-band. 

To get around the limitations of the currently available 
instruments, some have utilized the powerful technique 
of median stacking. By taking a population of objects 
selected at a different wavelength and combining cutouts 
extracted from a large but shallow survey, it is possible to 
probe far below the original survey threshold and deter- 
mine median trends in the data. For example, this tech- 
nique has been applied to the FIRST survey to study var- 
ious source popula tions with median fl ux densities in the 
10s - 100s of fiJy. iWhite et all (|2007f ) studied the radio 
properties of qu asars selected froni the Slo an Digital Sky 
Survey (SDSS; lAbazaiian et al.l I2OO90 . Ide Vries et all 
(|2007D applied the method to SPSS-selected active galac- 
tic nu clei (AGN) and star-forming galaxies. iHodge et al.l 
(HOOI) studied the radio emission from galaxies l acking 
strong optical emission lines, and Hodge et al.l ()2009f ) 
used median stacking to study low-luminosity AGN in 
luminous red galaxies (LRGs) out to redshifts of 0.7. 

This method is therefore a way of milking the most 
out of today's radio surveys, but, like everything, it too 
has its selection biases. In this case, the obvious bias is 
the limitation to samples selected at wavelengths other 
than the radio. To study certain faint radio populations 
without this bias, what is needed is a survey that is both 
fairly wide and fairly deep. Here we present a 1.4 GHz 
VLA survey of the SDSS Southern Equatorial Stripe, 
also known as "Stripe 82"0. Our survey is intermediate 
in both size and depth, covering over 90 deg^ to almost 
three times the depth of FIRST, or an rms of 52 /.tJy. 
This makes it the deepest 1.4 GHz survey to cover this 
much area, and conversely, the widest survey to reach 
this depth. 

On top of its unique position in depth versus width 
parameter space, this survey brings with it a significant 
amount of complementary data. SDSS Stripe 82 has ex- 
isting multi-wavelength coverage from the radio to the 
ultraviolet, the most noteworthy being SDSS imaging 
two magnitudes deeper than the typical SDSS observa- 
tions. Other coverage includes YJHK imaging by the 
UKIRT Infrared Deep Sky Survey (UKIDSS) as part of 
its Large Area Survey (LAS), and millimeter imaging by 
the Atacama Cosmology Telescope (ACT). The stripe 
also has multi- epoc/i coverage, with over 100 epochs of 
SDSS imaging and two previous 1.4 GHz epochs from 
the FIRST survey. Finally, our new VLA survey was 
done with 1.8" angular resolution, giving high- resolution 
radio imaging over the entire 90 deg^ area, as well as 
serving as a test project for a future FIRST-sized (or 

3 www.physics.drexel.edu/~gtr / vla/stripe82/ 



larger) high-resolution survey. 

All of these characteristics combine to make this sur- 
vey a powerful dataset for investigating the extragalac- 
tic radio sky; in particular, the AGN. Radio AGN have 
typically been divided into radio-loud and radio-quiet 
objects, with t he dividing line g enerally set near 10^^ 
W Hz-i, as in iFanaroff fc RilevI (|T97l . The cause for 
the dichotomy, and even its exi stence in the first place , 
has long been sub ject to debate (jKellermann et allll989l : 
llvezic et al.ll2002l ) . Recent studies seem to imply that the 
distribution is formally bimodal, though with a larger 
radio-intermediate population than previously thou ght 
(|Cirasuolo et al.l[2003t lLaMl[2003t IWhite et al.l[200l . 

There is also the question of the behavior of the di- 
chotomy as a function of redshift and luminosity. This is 
a tricky question to answer, as degeneracy between red- 
shift and luminosity inherent in flux limited surve ys may 
disguise a dependency. According to llvezic et al.l (|2b02i) , 
the fraction of objects classified as radio-loud or radio- 
quiet does not evolve wit h either red s hift o r luminosity. 
A more recent study by iJiang et all (|2007D found that 
the radio-loud fraction of quasars is a strong function of 
both redshift and luminosity, decreasing with increasing 
redshift and decreasing luminosity. 

Related to the evolution of the dichotomy is the evo- 
lution of the populations themselves. It has been known 
for some time that high-luminosity radio sources show 
strong cosmic evolution between z "^ and z ^ 2 
(iLongaii] [l966l : iDunlop &: Peacockl [l990l : iWillott et afl 
I2001D. The evolution is m ore uncertain beyond z ^ 
2 (IJarvisfcRawlinesI [20001) ■ and some studies suggest 
that it declines fe.g.. iWall et al.l [20051 ). The evolution 
of low-luminosity sources is also not well constrained. 
Their evolution seems to be weaker than that of pow- 
erful radio sources, and some studies imply very weak 
to almost no evolu tion at all (jC lcwlcv & Jarvis 120041: 
Sadler et ahl 120071: iDonoso et ai.i ,2009^ .Smolcic et al.l 



2009: IMcAlpine fc Jarvisll2010i) . 

Often coming hand-in-hand with measuring radio 
source evolution is the determination of the radio quasar 
luminosity function (QLF). Since it is now thought that 
AGN feedback ma y be crucial in the f ormation of mas- 
sive galaxies (e.s.. iHopkins et aLll2006[ ). the QLF is one 
of the key diagnostics in galaxy formation models. While 
the QLF has been deter mined self-consistentl y from the 
optical to the hard X-ray (|Hopkins et al.lf2"007l ). the radio 
QLF has only bee n determined for relatively bright (> 
250 mjy) sources (|Wall et al.ll2005f ). limiting the ability 
to constrain galaxy evolution models. 

The large sample of faint radio sources this survey will 
provide will help shed light on the issues discussed above 
and more. In this paper - the first on the new VLA 
coverage of the stripe - we describe the basics of the 
survey and the results of various tests of quality con- 
trol. The format of the sections is as follows: We be- 
gin in Section [5] with a discussion of the technical de- 
sign of the survey. In Section [3l we give an account of 
the observations, followed by an outline of our data re- 
duction steps (Section . We then present the source 
catalog (Section [5]), including a brief summary of the 
source extraction algorithm and definitions of each col- 
umn in the catalog. We include a thorough discussion 
of image characteristics (Section [T]), including achieved 
depth, astrometry, photometry, and morphology. Com- 
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Fig. 1. — Multiwavelength coverage of SDSS Stripe 82. Panel (a) shows the high resolution VLA coverage of this work, where the 'X's 
mark the individual pointing centers. Panel (b) shows the GALEX coverage. Panel (c) shows the initial coverage by ACT. ACT imaging 
covering the entire stripe should be public by the end of 2011. Also shown are the areas targeted by the WiggleZ Dark Energy Survey, and 
the coverage of the VIRMOS-VLT Deep Survey (VVDS), the UKIDSS Deep Extragalactic Survey (DXS), and the CFHT Legacy Survey 
Wide field 4. Panel (d) shows the DRl UKIDSS coverage, with the approximate DR4 coverage indicated by the dashed line. A random 
sampling of SDSS photometry is shown in panel (e), where the density increases toward the Galactic Plane. 



pleteness is discussed in Section [51 We conclude with a 
section discussing matching to the SDSS and the result- 
ing optical identifications (Section [9]). Where applicable, 
we assume a flat Friedmann-Lemaitre- Robertson- Walker 
(commonly FRW) cosmology of a spatially isotropic, ho- 
mogeneous universe with Ho = 70 km s~^ Mpc~^, J7a = 
0.7, and Om = 0.3. 

2. TECHNICAL DESIGN 

Stripe 82 is the SDSS term for the region of sky that 
runs along the Celestial Equator in the southern Galac- 
tic cap. It was imaged repeatedly by SDSS during the 
months when the northern Galactic cap was not observ- 
able. As a result. Stripe 82 has much deeper imaging 
than the rest of SDSS, reaching to g '^24.5 in this roughly 
300 square degree area. 

The field also has extensive multi- wavelength coverage. 
In the UV, data from the GAL EX All-Sky (AIS) and 
Medium Imaging Surveys (MIS) (|Martin et al.ll2005[ ) ex- 
ist over much of the fie ld, and the UKIRT I nfrared Deep 
Sky Survey fUKIDSS: iLawrence et"aIll2007D has targeted 
Stripe 82 as part of their Large Area Survey (LAS). 



The VISTA Hemisphere Survey ()Arnaboldi et al.|[2007[) . 
reaching Jab — 21.2, has started taking data and will 
eventually include Stripe 82. Millimeter imaging with the 
Atacama Cosmology Telescope (ACT; Kosowsky 200^ 
reaches a depth of 1 mJy and should be public by the 
end of 2011. 

Parts of Stripe 82 were als o previously covered in the 
radio by the FIRST survey (|Becker. White, fc Helfandl 
119951 hereafter BWH95), which has a typical rms noise 
of 0.15 mJy beam~^ and a resolution of 5". The large 
amount of multi-wavelength data made it an obvious 
choice for deeper radio follow-up. Our technical goal 
was therefore to produce a radio survey of Stripe 82 
to 3-5 times the depth of the FIRST survey. The data 
were taken primarily with the VLA's A-configuration to 
achieve the highest angular resolution possible. 

As the primary beam is the same size as that used for 
th e FIRST survey, we utilized the same pointing grid 
as IBWH95I . This grid keeps the right ascension spacing 
fixed within strips of constant declination, and adjusts 
the declination spacing as necessary so that the sensitiv- 
ity will be uniform once adjacent fields are coadded. In 
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this region {6 = 0°), that means pointings are separated 
by 3 min in RA and 13' in dechnation. 

Figure [T] shows the area covered with our high- 
resolution VLA observations, as well as some of the 
multi-wavelength coverage. The 'X's in the VLA map 
in Panel (a) show the individual pointing centers of this 
work. We avoided the region around RA = as infrared 
observations at wavelengths longer than '^bfiin would 
be affected by the zodiacal light band and contamina- 
tion by asteroids, which makes it a less useful multi- 
waveleng th region. Pane l (b) shows the GALEX cov- 
erage (Martin et al. 120051) . Panel (c) shows the coverage 
by ACT (iKosowskvl[200l . as well as the areas targeted 
by th e WiggleZ Dark Energy Survey (jPrinkwater et ahl 
l2010f ). Also shown by the small boxes is coverage by 
the V IRMOS-VLT Deep Survey fVVDS: EeFevre et al.l 
200l. the UKIDSS Deep Extragalactic Survey (DXS; 



Lawrence et"aLl I2007D . and the CFHT Le gacy Survey 



Wide fieldfl A random sampling of SDSS photometry 
is shown in panel (d), where the density increases toward 
the Galactic Plane. 

3. OBSERVATIONS 

The data were collected over two VLA cycles (2007- 
2008, 2008-2009). The majority of the observations were 
taken in the A-configuration, but we also obtained B- 
configuration coverage of the area to improve the sam- 
pling of the Fourier (U-V) plane and increase sensitivity 
to extended structure. In the first cycle, we were awarded 
scheduled A-configuration time (Program ID AR646), 
and the B-configuration time (AR659) was all collected 
dynamically through the queue. In the second cycle, we 
were awarded both the A and B-configuration time as 
dynamic time (AR685). We covered Area 1 (Figure [1]) in 
the A & B-configurations in 2007-2008, and Area 2 (in A 
& B-configurations) in 2008-2009. Area 1 is made up of 
275 pointings, and Area 2 has 374, coming to 649 fields, 
and 92 deg^ covered in total. 

Due to various issues with the data collection and 
the use of dynamic time for all but the 2007-2008 A- 
configuration, the total time on each field varies consider- 
ably. For the 2007-2008 scheduled A-configuration time 
(AR646), each field has anywhere from five to ten 284- 
second observations (time on source, or "dwell" time) 
devoted to it, or 1420-2840 seconds. For the 2007-2008 
B-configuration time (dynamic, AR659), we obtained at 
least one 284-second observation of each field, and as 
many as four 284-second observations (1136 seconds to- 
tal dwell) near 5 = 0° in Area 1. 

For the 2008-2009 cycle, we were awarded dynamic 
time in both the A and B-configurations. For the A- 
configuration, we divided our allotted 188 hours such 
that every field was observed for at least three 495-second 
intervals, for a total of 1485 seconds (dwell) per field. 
Some fields were affected by WIDAR testing, which uti- 
lized the easternmost antenna. Since good resolution 
is one of our goals for this project, we have also re- 
observed those fields during non- WIDAR time. For the 
B-configuration, we were not able to use all of our time 
due to the popularity of our LST range during this cy- 
cle. Some data had to be collected during the BnA- 
configuration. In the end, we succeeded in obtaining at 

* http:/ /www.cflit.hawaii.odu/Science/CFHLS/ 



least 220 seconds (dwell) on every field. 

To minimize bandwidth smearing, where sources away 
from the beam center are radially blurred due to the fi- 
nite receiving bandwidth, our observations utilize band- 
width synthesis. We observed in spectral line mode with 
7 frequency channels of 3.125 MHz each. Correlator lim- 
itations then set the bandwidths of the 2 IFs (for the 
two circular polarizations) at 25 MHz, and these were 
centered at 1465 and 1385 MHz. Our chosen integration 
time was 3 seconds. We used a number of phase calibra- 
tors spread out over the area in order to minimize slew 
time, and we typically observed a phase calibrator every 
30-60 minutes for 40-90 seconds. We tied the observa- 
tions to a flux scale using 3C48 as our flux calibrator, 
which was observed once during every observation block. 

4. DATA REDUCTION 

We reduced the data using the Astronomical Imaging 
Processing System (AIPS) software. As the data were 
taken during the transition from the VLA to the EVLA, 
we created a new Channel dataset from the line data 
in order to avoid closure errors. We flrst applied baseline 
corrections and ran UVFIX to correct post-MODCOMP 
uvw errors: some of our early data had u'w, v's, and 
w's written with the wrong sign, and there was also a 
problem with the online calculation of negative u's, v's, 
and w's. Then we flagged the data by hand with CLIPM 
and UVFLG to remove any problems due to interfer- 
ence, cross-talk between antennas, correlator failures, or 
bad antennas. All the data went through four rounds 
of self-calibration (two phase-only, and two phase and 
amplitude) . 

We imaged the data with IMAGR using multiple facets 
and three-dimensional techniques. We used 265 facets 
with 1024 0.3" pixels per side. The flelds used for clean- 
ing were thus very large, and we did not find it necessary 
to take into account sources outside the field. The IFs 
were imaged together in all fields except a small num- 
ber of problem fields where the resultant rms was sig- 
nificantly higher than the survey median. To keep as 
uniform as possible, we cleaned each field down to 0.2 
mjy. There is also a trade-off between sensitivity and 
resolution. Due to the declination of this stripe, the VLA 
beam is significantly extended in the N-S direction from 
the nominal A-configuration resolution of 1.4". To main- 
tain a circular beam while using natural weighting, we 
found that the resolution degraded to 2.5". In the final 
round of imaging, we therefore use uniform weighting 
with robust=0 and a circular beam of 1.8". The facets 
were stitched back together with FLATN using a cutoff 
radius of 22.3', and the final images are 4650 x 4650 0.6" 
pixels. These images, which were made from individual 
pointings, are the so-called "grid" images. The typical 
rms for the grid images is 52 fiJy (see Section [7T|) . 

To produce the final maps, we co-added the individual 
grid images using an adaptation of the pipeline produced 
for the FIRST survey by BWH95. This step increases 
the S/N and improves sensitivity away from the beam 
center. The pipeline searches for fields in the vicinity of 
the map center, truncates each where the primary beam 
sensitivity falls below 18%, weights by the square of the 
primary beam response, and sums the fields. Due to 
overlapping coverage, we produce co-added maps for only 
every other declination strip. The resultant maps are 
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4650 X 4650 pixels, or 46.5' x 46.5' in extent, with 0.6" 
pixels. 

5. THE CATALOG 
5.1. Source Extraction 

To create the catalog, we used HAPPY, the source- 
finding algorithm which was developed initially for use by 
the FIRST survey I Whiteltall (fTOOl . HAPPY searches 
an image for pixels exceeding a user-defined threshold, 
and it places rectangles around contiguous regions of such 
pixels to form "islands" . Each island is then fit with up to 
four two-dimensional Guassian components, which must 
pass several criteria for acceptance. The end result is a 
list of elliptical Gaussia n components . Mor e details on 
HAPPY can be found in I White et all ()1997| ). 

HAPPY requires a noise threshold to be specified, 
which we compute as a function of position in a coverage 
map. To create the map, we combine the rms HAPPY 
computes by looking at a histogram of pixel values for 
each grid image with a model for the elevated noise 
around bright sources, with the noise being much higher 
in the north-south direction. This additional step is nec- 
essary in order to ensure the strong sidelobes present 
around many of the bright sources in this survey do not 
end up in the source catalog. We then add the contribut- 
ing noise maps at each position, weighting them in the 
same way we weight the grid images that make up the 
coadded images. We set HAPPY's noise threshold as five 
times the rms found in this coverage map at the position 
of interest. To view the coverage map for Stripe 82, see 
Section O 

Since the coadded images have some overlap, a catalog 
produced in this way will have duplicate sources. Iden- 
tifying the duplicates is generally straightforward, since 
the images that contribute at any sky position are ex- 
actly the same in the coadded images that cover that 
position. The only difference is that the pixel grid cen- 
ters will be slightly different. For the vast majority of 
sources, the fitted parameters are almost identical, and 
they agree much more closely than the actual uncertainty 
in the parameters. 

There may be cases, however, where a source gets 
decomposed into multiple components, and the source 
fits are therefore different in different fields. To identify 
groups of duplicate sources and select which sources to 
retain from the matched groups, we utilize the following 
algorithm: We first search for matches within 1.7" of each 
source, where 1.7" is chosen because it is approximately 
the size of the PSF FWHM. Groups may be chained to- 
gether by links of <1.7", so the total group size may be 
larger than this if there are multiple links. Sources that 
are separated by more than 1.7" are resolved into sep- 
arate components. We then find all coadded fields that 
overlap the mean RA/Dec position for the group, and 
we identify the field where the group is closest to the 
center as the best field. We then keep only the sources 
that come from that best field, discarding duplicates from 
fields other than the best one. This is differ ent than the 
appro ach used by FIRST and outlined in iWhite et al.l 
(|1997D . and it is applied to all sources, not just those 
with duplicates. The final catalog should therefore only 
contain sources from fields that are identified as the best 
field at that source position. Applying this process to the 



Stripe 82 data, we produce a final catalog of 17969 radio 
sources. An example page from the catalog is shown in 
Table [B 

5.2. Column Definitions 

The columns are well-defined on the FIRST survey 
websit^B, but we give a brief summary here. 

Columns (1) and (2) - Right ascension and declina- 
tion of the source in J2000. The positional errors are 
a function of source brightness, size, and noise in the 
map. They are best found using a simple rule-of-thumb 
approach, as the HAPPY-derived errors tend to be un- 
derestimated. An empirical equation for the accuracy at 
90% confidence is /Sizex {1/ SNR+1/2Q), where /Size is 
the fitted major o r minor a xis size, and SNR is the sig- 
nal to noise ratio (t White et al. ,1997). Systematic errors 
are smaller than 0.05". 

Column (3) - P(S) indicates the probability that the 
source is spurious, most likely because it is a sidelobe of 
a nearby bright source. Low values of P(S) mean the 
source is unlikely to be spurious. The probabilities are 
computed using an algorithm based on multiple voting 
oblique decision tree classifiers, which were trained on 
deep VLA fields. The algorithm was developed for the 
FIRST survey and will be described in more detail in a 
future paper on the final FIRST catalog. Note that the 
algorithm is optimized for the FIRST survey, whereas the 
RMS computation for this catalog has changed signifi- 
cantly. The values of P(S) are therefore not very reliable 
for this catalog. 

Column (4) - Fpeak is the peak fiux density in mJy 
beam~^ derived by fitting an elliptical Gaussian model 
to the source. The uncertainty is given by the rms noise 
at the source position. 

Column (5) - Fint is the integrated flux density in mJy 
derived from the elliptical Gaussian model fit. The un- 
certainty in Fint can be considerably greater than that 
of Fpeak depending on source size and morphology. An 
expression to est i mate the uncertainty can be found in 
iSchinnerer et al.l (|2004[ ). For point sources, the relative 
uncertainty (c///) reduces to: 

^ = ^2.5^+0.01^ 

Column (6) - RMS is the local noise estimate at the 
source position in mJy. The RMS is calculated by com- 
bining the measured noise from all grid images contribut- 
ing to the coadded map at the source position. 

Columns (7), (8) and (9) - Maj, Min, and PA give the 
major and minor axes (FWHM in arcsec) and position 
angle (degrees east of north) derived from the ellipti- 
cal Gaussian model for the source. Maj and Min have 
had the elliptical Gaussian point-spread function decon- 
volved. Noise can cause the fitted values of major and 
minor axis prior to deconvolution to be smaller than the 
beam, and the deconvolved size is given as zero in that 
case. The uncertainties in the deconvolved sizes depend 
on both brightness and size. 

Columns (10), (11) and (12) - fMaj, fMin, and fPA 
give the major and minor axes (FWHM in arcsec) and 
position angle (degrees east of north) derived from the 

^ http://sundog.stsci.edu 
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TABLE 1 

Excerpt from the VLA Stripe 82 Survey Catalog 



SDSS Matches 

RA Dec P(S) Fpeak Tint RMS Maj Min PA £Maj fMin fPA Field # Sep i CI 

(J2000) (J2000) (mJy) (mjy) (mjy) (") (") (deg) (") (") (deg) (") (mag) 



22 26 25.132 


+00 


25 09.26 


0.317 


0.52 


0.29 


0.092 


0.00 


0.00 


22 20 19.710 


+00 


25 08.57 


0.276 


0.53 


0.45 


0.085 


0.00 


0.00 


02 06 42.812 


+00 


25 08.36 


0.200 


0.97 


1.32 


0.063 


1.14 


1.01 


22 11 36.231 


+00 


25 07.48 
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''NOTES. See Section 5.2 for column definitions. 

elliptical Gaussian model for the source. These are the 
fitted sizes before deconvolution of the 1.8" circular clean 
beam. 

Column (13) - The field name is the name of the 
coadded image containing the source. Note that the 
field name encodes the center of the field; field hh- 
mmm+ddmmm is centered at RA = hh mm.m, Dec = 
+dd mm.m. The letter appended to the field name in- 
dicates the last catalog release in which the image was 
modified. 

Columns (14), (15), (16), and (17) give information 
on optical counterparts from the SDSS DR6, which were 
found using the CasJobs web interface. Column (14) 
gives the number of matches within a radius of 8". For 
the closest match. Columns (15) - (17) give (respec- 
tively) the separation in arcsec, magnitude in the SDSS 
i-band, and morphological class, where s=stellar and 
g=nonstellar /galaxy. 

6. IMAGE ANALYSIS 

This section introduces two effects, bandwidth smear- 
ing and resolution bias, which will influence the derived 
catalog properties discussed in the next section (Section 
[7|) in complex ways. For further discussion of these effects 
and possible correction schemes, the reader is referred to 
(for c.g^ Huvnh et al. (2005.), .Prandoni et al.. (,2001.) . 
and iBondi et all (2008D . 

6.1. Bandwidth Smearing 

The effect of bandwidth smearing is to radially-smear 
the source flux, thus causing objects to appear more ex- 
tended and the peak flux to be underestimated. To cal- 
culate the effect of bandwidth smearing on the beam size, 
we plot measured major axis size versus distance from the 
field center for all the grid image sources with Speak > 5 
mJy beam^^ (Figure [2]). We use only bright sources to 
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Fig. 2. — Fitted major axis size versus distance from field center 
for all grid image sources with Speafe > 5mjy. 

ensure the results are not affected by spurious sources at 
fainter fiux densities. Some sources are resolved, which 
produces a scattering of points toward high values of the 
major axis, but there is a clear lower bound that is de- 
fined by point sources. We fit a function to this ridge 
(red line) and used it to compute the weighted, smeared 
beam minor axis versus major axis for all the catalog 
sources. We do this by calculating the moments (x^, y^, 
xy) for the point spread function in each grid image in- 
cluding the effect of bandwidth smearing, which stretches 
the sources in the radial direction. Depending on where 
the source falls in the pattern of grid images, there may 
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Fig. 3. — The measured weighted, smeared beam minor and ma- 
jor axes for all the catalog sources. Without bandwidth smearing, 
all sources would be at the position (1.8", 1.8"). 
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Fig. 4. — The histogram of major and minor axis sizes for the 
points in Figure[3] The median major axis is 2.13", and the median 
minor axis is 1.95". 

be several contributing images with stretching oriented 
in different directions. We calculate the weighted sum 
of all these moments (in the same way we weighted the 
grid images to produce the coadded image), and convert 
back to a Gaussian model with the major axis, minor 
axis, and position angle computed from the moments. 

Figure [3] shows the weighted, smeared beam minor axis 
versus major axis for all the catalog sources. Without 
the effect of bandwidth smearing, all points would be at 
the position (1.8, 1.8). Due to bandwidth smearing, the 
bulk of the points have slightly large major and minor 
axes. (Note that since the beam may be smeared in dif- 
ferent directions in the different grid images, the minor 
axis is also typically larger than 1.8"). Complex struc- 
ture is visible in the Figure, which traces the geometry 
of the overlapping regions between fields. There is also 
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Fig. 5. — The solid histogram shows the distribution of values 
of peak flux density for FIRST sources that have no match in 
our new Stripe 82 radio catalog. The dotted histogram shows the 
distribution of values of peak flux density for all FIRST sources in 
the area, including those with matched sources in the new Stripe 
82 catalog. Both histograms have been normalized to a peak value 
of 1.0. 

an obvious scattering of points which extend out to even 
larger beam sizes (major axis '~3"). These points repre- 
sent sources at the edge of our coverage, where there are 
no observations close to field-center, and all observations 
are therefore far off axis. 

Figure HI shows the histogram of major and minor axis 
sizes for the distribution of points in Figure [31 The me- 
dian major and minor axis sizes are 2.13" and 1.95", re- 
spectively, and the means are essentially identical to the 
medians. There is a well-defined peak at a major axis of 
2.2". There is also a peak at exactly 1.8" corresponding 
to sources at the centers of grid images where there are 
no overlapping grid images, and therefore no bandwidth 
smearing. The current catalog version does not include 
a correction for the effect of bandwidth smearing. One 
possible fix would be to degrade everything to the same 
resolution by smearing the beams with a kernel of the 
appropriate size prior to coaddition of the grid images. 
It is important to keep in mind that bandwidth smear- 
ing will affect all of the survey characteristics discussed 
in Section [71 

6.2. Resolution Bias 

Resolution bias is due to the fact that resolved sources 
with integrated fluxes over the catalog cut will not make 
it into the catalog if their peak fiux is below the cut. 
As the Stripe 82 survey was conducted at high resolu- 
tion, we therefore expect to see a large effect from the 
resolution biasl^. One simple way to estimate the effect 
of the resolution bias is to compare our catalog to the 
lower-resolution FIRST survey. By looking at how many 
FIRST sources are undetected in the Stripe 82 catalog, 
we can get an estimate of the number of sources that are 
being "resolved out" between the FIRST catalog and the 
Stripe 82 catalog. 

^ One might ask why the catalog is not instead constructed so 
that integrated flux is the limiting criterion; however, as interfer- 
ometers will always "resolve out" flux from extended structures 
with angular sizes larger than the largest available resolution ele- 
ment of the array, a catalog constructed in this fashion would suffer 
from its own incompleteness. 
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Fig. 6. — The solid histogram shows the distribution of decon- 
volved major axis sizes for FIRST sources that have no match 
in our new Stripe 82 catalog. The dotted histogram shows the 
distribution of sizes for all FIRST sources in the area, including 
those with matched sources in the new Stripe 82 catalog. Both 
histograms have been normalized to a peak value of 1.0. Nega- 
tive values of deconvolved major axis indicate sources that have 
a fitted major axis size smaller than the beam and likely are due 
to sidelobes, which often appear elongated in one direction and 
narrow. 

Of the 11008 FIRST sources in the Stripe 82 area, 8635 
have one or more match within 5" in the Stripe 82 cat- 
alog, leaving 2373 FIRST sources without a Stripe 82 
counterpart (22%). A (normalized) histogram showing 
the distribution of values of peak flux density for these 
FIRST sources is shown in Figure [5] (solid line). For 
comparison, in the same figure we also show the distri- 
bution of peak flux density values for all of the FIRST 
sources (dotted line). The distribution of unmatched 
FIRST sources is clearly skewed toward low values of 
peak flux density. This result is consistent with what we 
would expect if resolution effects are at play, since fainter 
sources are more easily resolved out. 

To determine whether the unmatched FIRST sources 
tend to be extended, we plot the normalized distribution 
of deconvolved major axis sizes compared to that of all 
FIRST sources in Figure [51 The unmatched sources ap- 
pear to be slightly skewed toward larger sizes. As a pre- 
caution, we repeated the matching between the FIRST 
and new Stripe 82 catalogs, but using a matching radius 
equal to one-half of the fitted major axis size instead. 
The matching radius thus scales with source size, ensur- 
ing that we are not missing a valid population of resolved 
but legitimate sources by not searching far enough from 
the FIRST source center. This match still resulted in 
^21% of the FIRST sources lacking a Stripe 82 counter- 
part. We conclude that, as expected, extended sources 
are more likely to be resolved out. The resolution bias 
will affect the photometry, completeness, and morphol- 
ogy of the catalog, which will be discussed in the follow- 
ing section. 

7. SURVEY CHARACTERISTICS 
7.1. Depth 

Figure [7] shows a histogram of the rms noise levels for 
the 649 grid images that make up our co-added maps, as 
well as the cumulative fraction of fields as a function of 
rms noise. The median rms is 52 /zJy beam~^, and 97% 



of the fields have rms noise levels < 60 /i Jy beam ^ . We 
therefore achieve our target rms noise of three times the 
depth of FIRST, or ~50 /iJy beam'^ We note that 
there are a few images where the noise is determined by 
the dynamic range, and they are compromised by bright 
sources. We also note that in the immediate vicinity of 
bright sources, the noise is significantly worse than the 
total image rms. This will be discussed below. 

A coverage map, showing the rms sensitivity achieved 
in the co- added images, is shown in Figure [51 The rms 
sensitivity is shown for the entire region surveyed with 
a resolution of 1' x 1', which is only for presentation 
purposes, as the catalog is created by computing the rms 
at the exact position of each source. Note that the x and 
y axes are not to scale, as the surveyed strip is very thin 
in Declination. To make the map, we created a matrix 
of the desired resolution, and for each position in the 
matrix, we determine which grid images contributed to 
the co-added map at that position. We then took the 
rms values of the grid images and convolved them with 
a model for the noise around bright sources. This was 
necessary in order to exclude sidelobes from the final 
source catalog. 

Finally, we combined these values with the same 
weighting used to co-add the grid images. The resul- 
tant rms sensitivity values per 1 arcmin^ region range 
from 35 /iJy bearn"^ to 4.4 mJy beam~^. The high end 
is due to areas in the immediate vicinity of very bright 
sources. The sensitivity map is used to establish the 5a 
threshold for catalog sources. 

7.2. Astrometry 

To assess the astrometric accuracy of the final maps, 
we tested the robustness of our catalogued source posi- 
tions both internally and against external catalogs. To 
do the internal test, we utilize the fact that there is a 
fairly significant amount of overlap between adjacent grid 
images. A single source can therefore appear in multiple 
grid images, and calculating the positional offset between 
observations for each multiply-observed source is a good 
test of astrometry for weak sources. 

Figure |9] shows the positional offsets for multiply- 
observed sources with S > 1 mJy. Each point is plotted 
twice - the second time with the sign reversed - since the 
order is arbitrary. The ellipse shows two times the rms 
scatter in right ascension and declination. Note that we 
tried restricting the matching to point sources, but we 
found that the rms does not improve by putting a limit 
on deconvolved major axis. We therefore allow all sources 
in the matching. This plot shows 4598 pairs of sources 
and has an rms of 0.19" in both right ascension and dec- 
lination. For comparison, the same analysis for FIRST 
yields rms values of 0.49" and 0.46" in right ascension 
and declination ()BWH95D . This survey thus achieves 
an astrometric accuracy approximately 2.5 times better 
than the FIRST survey. Although the resolution of the 
A-configuration is three times better than FIRST, we be- 
lieve that our astrometry may be limited by bandwidth 
smearing, which was discussed in Section [6.11 

For the external checks, we compared against the 
FIRST survey and the SDSS DR7 Quasar Catalog 
([Schneider et al.[[2010D . (Note that we would have liked 
to compare against the VLA calibrator catalog, but there 
are only seven VLA calibrators in the area covered by 
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Fig. 7. — Left: a histogram showing the rms noise levels for the 649 grid images that make up our co-added maps. Right: the cumulative 
fraction of fields as a function of rms noise. 
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Fig. 8. — Sensitivity map showing the rms noise achieved across the surveyed area. Note that the x and y axes are not to scale. 



these observations — not enough for a statistically signif- 
icant comparison.) Figure ITOl shows the comparison with 
the FIRST survey. As fainter sources will have inherently 
worse astrometric accuracy, we limit the comparison to 
sources that are both point sources and have flux densi- 
ties greater than 10 mJy in the FIRST catalog. The 672 
source pairs represented in Figure [10] have rms values 
of 0.13" and 0.18" in SRA and SDec, respectively, with 
mean offsets of -0.006" ± 0.005", and -0.038" ± 0.007", 
one-fifth the rms accuracy of overlapping sources. The 
ellipse again shows two times the rms scatter in right 
ascension and declination. These values of rms scatter 
are no better than the FIRST survey, which had lower 
resolution; however, this merely reflects the fact that we 
are comparing against FIRST and therefore are limited 



by its accuracy. There does appear to be a slight offset 
in declination, in the sense that the Stripe 82 positions 
are biased slightly south of the FIRST positions. It is 
unclear what is causing the plume of points outside the 
ellipse at -0.15, -0.4, though they seem to occur in lines 
on the sky. This may be due to a small systematic cali- 
bration error in the positions. 

Figure [TT] shows the comparison between the Stripe 
82 catalog and the SDSS DR7 Quasar Catalog 
([Schneider et al.l[2010H for sources brighter than 5 mJy 
in the FIRST catalog. Such a comparison can be used to 
tie together the radio and optical reference frames. The 
75 source pairs have rms values of 0.11" and 0.14" in 
SRA and 5Dec, respectively, with mean offsets of 0.02" 
± 0.01", and 0.01" ± 0.02". Here there appears to be 
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Fig. 9. — The positional offsets for multiply-observed sources 
with S > 1 mjy. The ellipse shows two times the rms scatter in right 
ascension and declination. Our high-resolution survey achieves an 
astrometric accuracy approximately 2.5 times that of FIRST, with 
an rms scatter of 0.19" in right ascension and declination. 



Fig. 11. — A comparison of Stripe 82 source positions with the 
SDSS DR7 Quasar Catalog (Schneider et al. 2010). The compari- 
son was limited to sources brighter than 5 mJy. The 75 source pairs 
have rms values of 0.11" and 0.14" in 5RA and (5Dec, respectively, 
with mean offsets of 0.02 ± 0.01, and 0.01 ± 0.02. The rms scatter 
is very good, and there appears to be no offset in right ascension 
or declination. 
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Fig. 10. — A comparison of our source positions with FIRST 
sources brighter than 10 mJy. The ellipse shows two times the rms 
scatter in right ascension and declination. The values of rms scatter 
are consistent with that expected due to the accuracy of the FIRST 
survey. There does appear to be a slight offset in declination, in 
the sense that the Stripe 82 positions are biased slightly south of 
the FIRST positions. 

no significant offset in eitlier riglit ascension or declina- 
tion. We conclude that the astrometric accuracy of the 
Stripe 82 catalog is excellent, as confirmed by both in- 
ternal checks and comparisons to independent catalogs. 

7.3. Photometry 

We next check the photometric accuracy of the Stripe 
82 catalog. There are multiple sources of photometric 
error, including resolution bias, CLEAN bias, the model 
of the primary beam, and bandwidth smearing. We first 
compare our integrated flux densities with those from 
the FIRST survey for sources in both catalogs. This is 
shown in Figure [12] We restricted the matching to point 



sources with FIRST Smt > 1 mJy and with only one 
source per 'island', which is a concept used by the source- 
finding algorithm 'HAPPY' and is basically a contiguous 
set of pixels above a user-defined threshold (jWhite et al.l 
[1993). This Fi gure shows that there is a bias toward 
higher Stripe 82 values at low flux densities despite the 
fact that the Stripe 82 observations may resolve out ex- 
tended fiux more so than FIRST. This offset is likely due 
to the Malmquist bias. For sources that are unresolved 
in the Stripe 82 data, there would not be any bias if the 
surveys were equally sensitive. For sources that are re- 
solved in the higher-resolution Stripe 82 data, however, 
the Stripe 82 data is not as sensitive as the FIRST sur- 
vey. A source that is small compared with the FIRST 
beam (< 3-4" FWHM) but is resolved in the Stripe 82 
data will have a lower peak flux density in the Stripe 
82 data, which can put it below the Stripe 82 detection 
threshold. It will appear in the Stripe 82 catalog only if 
a noise fluctuation causes the apparent peak flux to rise 
above 1 mJy, in which case it will also appear brighter 
than the matched FIRST source. 

To test whether this is correct, we have plotted the dif- 
ference between the Stripe 82 and FIRST integrated flux 
densities as a function of Stripe 82 deconvolved major 
axis. This plot is shown in Figure 1131 We have limited 
the plot to weak (S < 3 mJy) sources, where the effect 
should be more prominent. We see that for unresolved 
sources, with a deconvolved major axis around 0.0", the 
distribution of points is fairly symmetric about zero. As 
the source size increases, however, the distribution shifts 
up toward higher Stripe 82 fluxes as the fainter sources 
disappear. This confirms that the Malmquist bias is at 
play. 

Another way to test for errors in the photometry is to 
again utilize multiply-observed sources. The idea is that 
if photometric errors occur in a certain field, then after 
the primary beam correction, the error will be a function 
of distance from field center. We show fiux densities for 
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Fig. 12. — A comparison of our integrated flux densities with 
those from FIRST for sources in both catalogs. This Figure shows 
that there is a bias toward higher Stripe 82 values at low flux densi- 
ties, which is likely caus ed by the Malmquist bias — see Section [7.3l 
Also refer to Figure [T3l which examines the difference between flux 
densities as a function of source size. 

pairs of multiply-observed sources in Figure fHl The "pri- 
mary field" refers to the field where the source is closest 
to the center. To make the fluxes comparable, we have 
corrected the peak flux of the primary field by the ratio 
of the primary beam corrections for the two fields. This 
is thus also a test of the primary beam model. Figure [14] 
shows a very tight correlation over almost three orders 
of magnitude, indicating that the relative photometric 
errors are small. 

Outliers in such a plot could also indicate variability, as 
the observations of the fields were taken anywhere from 
minutes to weeks apart. The variability would have to 
be strong, since each field is the combination of several 
individual observations and this would tend to damp out 
short term variability. There appear to be several promi- 
nent outliers in this Figure, all fainter than 1 mJy in the 
primary field. Upon closer inspection, however, these 
sources all show up in fields that are unusually noisy or 
have wave-like patterns due to incompletely-removed in- 
terference. If we throw out those sources, the width of 
the remaining correlation still implies we are only sensi- 
tive to variability of a factor of several on the faint end. 
We will conduct a search for variability in a follow-up 
paper (Hodge et al., in prep.). 

All VLA snapshot images suffer from a "CLEAN 
bias", a photometric defect where the fiux den- 
sity of sources is sys t emat ic ally underestimate d 
(IBecker. White, fc Helfandl [l99l ICondon et all [19981) . 
While our maps were not created from snapshot images, 
which should preclude a CLEAN bias, we searched for 
a systematic photometric bias by planting fake sources 
in the calibrated u-v data and repeating the exact same 
processing steps on the data to recover the sources. 
The planted sources were inserted with the AIPS task 
'UVMOD' and consisted of point sources with a range 



of peak flux densities and distances from field center. 
We randomly selected thirteen fields in which to insert 
the fake sources, and these fields had rms noise values 
between 45 and 73 /xJy beam^^. After recovering 
the planted sources with HAPPY, we compared the 
measured flux density of each source with its true flrrx 
density. 

In Figure [15] (left), we show the difference between 
measured and true flux density as a function of planted 
source strength. For each value of source flux density, we 
inserted twenty sources in a circle of diameter 12' in the 
u-v data. The data points at the top of the plot repre- 
sent the bias calculated using peak flux, while the data 
points at the bottom represent the bias calculated using 
integrated fiux. The bias is clearly small (roughly 35 /iJy 
beam~^ for peak flux density and sources 12' from the 
beam center) and it does not appear to depend strongly 
on source strength. 

To investigate the radial dependence of the bias seen in 
our Stripe 82 images, we inserted fake 1 mJy sources in 
the u-v data at four different distances from fleld center. 
The result, shown in Figure [T5] (right), is that the bias 
does not seem to depend on distance from field center. 
At all distances, it maintains a value of roughly 35 /iJy 
beam~^ as derived from values of peak fiux. This is dif- 
ferent than the CLEAN bias, where the functional form 
describing this effect happened to be such that, once the 
grid images were corrected for the shape of the primary 
beam prio r to being co-added, the radial de pendence dis- 
appeared ([Becker. White, fc Helfandl ll 9951) . This com- 
plicates matters somewhat, since the correction factor of 
a particular source would be a function of its position in 
each grid image that went into the co-added map. Be- 
cause of this, and since the bias is relatively small, we do 
not correct for it in the resultant Stripe 82 catalog. The 
result is that the bias for the grid sources will be 35 /^Jy 
beam~^ at field center and increase monotonically out- 
ward. Grid sources far from the field center will therefore 
see larger values of the bias. Note, however, that when 
the grid images are co-added, these sources will also con- 
tribute less to the co-added source flux density due to 
weighting by the primary beam. 

The largest effect on the photometry is therefore likely 
due to bandwidth smearing. The effect of bandwidth 
smearing is to reduce the peak flux density of the sources 
while keeping the integrated fiux density the same. Fig- 
ure [2] (Section 16. ip shows the measured major axis size 
versus distance from the fleld center for all the grid image 
sources with Speak > 5 mJy beam"-'^. This Figure shows 
that the smearing is relatively small out to 10'. Using the 
dimensionless parameter (3 = iiL—i^ — as an indicator 
(see [Bridle &: Swablll999f ). source peak flux densities are 
expected to be reduced by <10% out to that radius. This 
is an additional source of uncertainty that is not reflected 
in the generic equation for the uncertainty of catalog flux 
densities. 

7.4. Morphology 

One of the main differences between this survey and 
the FIRST survey is the approximately threefold increase 
in angular resolution. In addition to the enhanced detail 
provided on all sources, some sources are now resolved 
into multiple sources. To estimate the effect of this in- 
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Fig. 13. — The difference between the Stripe 82 and FIRST integrated flux densities as a function of Stripe 82 deconvolved major axis. As 
the source size increases, the distribution shifts toward higher Stripe 82 fluxes as the fainter sources disappear, implying that the Malmquist 
bias is at play. Note that these are weak sources. 
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Fig. 14. — A comparison of flux densities for pairs of multiply- 
observed sources. The "primary fleld" refers to the fleld where the 
source is closest to the center, and it includes a correction factor to 
account for the radially-dependent sensitivity of the primary beam. 

crease, we measured how many FIRST sources were re- 
solved into multiple components in our Stripe 82 radio 
catalog by matching the catalogs with a matching ra- 
dius of 5". We found that 1302 out of the 9328 de- 
tected FIRST sources were resolved into multiple com- 
ponents, or roughly 14%. Some examples of such sources 
are shown in Figure [T51 as 1' x 1' postage stamps. The 
images in the top row are the images from this survey, 
while the bottom row shows the corresponding cutouts 
from the FIRST survey. These sources all correspond to 
known SDSS quasars. 



The cumulative distribution function of major axis 
FWHM values for all sources in the catalog is shown in 
Figure [iTl The inset shows the function for small values 
of the major axis {"ip < 4"). The high fractions achieved 
for these small values illustrate that over 17,300 sources, 
or 96% of the catalog, are larger than the CLEAN beam 
size of 1.8". This does not represent the fraction of re- 
solved sources, as we established in Section that the 
effective beam size is actually larger due to bandwidth 
smearing. We use 2.3" as an estimate of the major axis 
value below which the vast majority of smeared beams 
lie (see Figure H]). Sources at the edge of the cover- 
age will have smeared beams larger than this value, but 
many sources will also have smeared beams somewhat 
smaller. Using this value, we find that ^60% of the cat- 
alog sources have been resolved. For comparison's sake, 
only ~20% of the sources in the FIRST catalog are larger 
than the FIRST CLEAN beam of 5.4". Our value is com- 
parable to other 1.4 GHz surveys of similar resolution, 
which probe much fainter flux densities (rms of a few 
/zJy) and find that '^70% of their sources are resolved 
when they adopt a classification method utilizing the 
resul t s of JMFIT (Owen fc Morrison i2008t IMiller et all 
[200llSchinnerer et al.ll20M ! 

8. COMPLETENESS 

The are several factors that affect the completeness of 
a source catalog. The most obvious obstacle to catalog 
completeness in our case is the fact that the rms sensi- 
tivity varies over the area covered, and sometimes quite 
dramatically (Figure [8] and Section I7.1|) . In addition 
to varying rms sensitivity, there are several additional 
factors that also influence the completeness, including 
CLEAN bias, bandwidth smearing, and the resolution 
bias. Here, we will compare the Stripe 82 catalog to ex- 
isting catalogs in order to estimate the completeness of 
the survey. 

In Section 16. 2[ we matched the Stripe 82 catalog to 
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Fig. 15. — Top: The difference between measured and true flux density as a function of artiflcial source strength. For each value of source 
flux density, we inserted twenty sources in a circle of diameter 12'. The data points at the top of the plot represent the bias calculated 
using peak flux, while the data points at the bottom represent the bias calculated using integrated flux. Bottom: The difference between 
measured and true flux density as a function of radial distance to the fleld center. 
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Fig. 16. — The top row shows 1' x 1' images resulting from this survey of several known SDSS quasars. The bottom row shows the cor- 
responding region of sky as seen by the FIRST survey. The columns show (1) J011129-I-003431, (2) J015243-I-002039, (3) J221308-I-003655, 
and (4) J222744-I-005730. 



the FIRST survey catalog and found that 22% of FIRST 
sources have no counterpart in the Stripe 82 catalog. As 
extended sources are more likely to be resolved out when 
they are faint, the completeness increases with decreas- 
ing source size and increasing source strength. To further 
understand the FIRST sources missing from the Stripe 82 
catalog, we compared their optical match rate with the 
match rate of FIRST sources that were recovered in the 
Stripe 82 catalog. For the optical data we used the deep 
SDSS imaging of Stripe 82, which is available through 
CasJobfl Cross-matching the catalogs with a matching 
radius of 5" reveals that FIRST sources with Stripe 82 
radio counterpart (s) have an optical match rate of 72%, 
whereas FIRST sources without Stripe 82 radio coun- 
terparts have an optical match rate of only 46%. The 



Specify Context = "Stripe 82" and runs 106 or 206 to retrieve 
the deeper co-added photometry. 



background rate is unimportant, as it would be the same 
for both scenarios. This drop in match rate therefore 
indicates that some of the undetected FIRST "sources" 
may not have been sources at all, but noise peaks and im- 
age artifacts. Based on the percentage drop, it is possible 
that as many as 26% of the unmatched FIRST sources 
may have been false detections. We checked the sky dis- 
tribution of unmatched sources, and some of them do 
group into narrow lines, a pattern indicative of sidelobes. 
However, after examining some of the unmatched sources 
by eye, we conclude that the percentage of false sources 
is probably not quite this high. While some of the un- 
matched sources are indeed sidelobes, the vast majority 
are undetected for other reasons, typically because they 
have been completely resolved out, or are visible but just 
below the Stripe 82 detection threshold. Most are faint 
sources right at the FIRST detection threshold of 1 mJy. 
We next looked to the population of previously-known 
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Fig. 17. — The cumulative distribution function of major axis 
FWHM values for all Stripe 82 sources. The inset highlights small 
values of the major axis up to ip = 4". 

radio quasars to see how many were recovered in the 
Stripe 82 catalog. This population should provide a sam- 
ple of sources where (1) the objects are not expected to 
be significantly extended, and (2) the objects are con- 
firmed as real sources. The re are 3885 quasars fr om the 
SDSS DR7 Quasar Catalog (jSchneider et al.ll2010D in the 
area covered by our observations. The FIRST catalog has 
sources matching 229 out of the 3885 within a matching 
radius of 5". Of those 229, the Stripe 82 radio catalog re- 
covers 223 of them. Of the 6 missing, all are visible by eye 
in the cutout server at the position of the FIRST sources. 
Five of them appear to be barely resolved out, with faint 
peak flux densities just below the catalog threshold at 
their position. The remaining source is brighter, with 
a peak flux density approaching 1 mJy, but it has been 
excluded from the catalog due to a bright source in the 
vicinity, requiring the rms threshold to be automatically 
raised to avoid adding spurious sidelobes to the source 
catalog. The new catalog is therefore highly complete for 
the radio quasar population, recovering over 97% of the 
FIRST-detected quasars. 

The discussion so far has been limited to a compari- 
son with the FIRST survey, which has a detection limit 
of only 1 mJy. To estimate the Stripe 82 survey com- 
pleteness below 1 mJy, we will compare our 1.4 GHz 
sour ce counts with th ose from the VLA-COSMOS sur- 
vey (jBondi et al.ll2008[ ) , which was also conducted at high 
resolution (1.5"), but achieves a much lower rms sensi- 
tivity of 11 iiJy. We calculated the Stripe 82 differential 
counts by dividing the number of sources in each flux 
density bin by the survey area of that bin and the bin 
width. We normalized the differential counts using the 
historical value of S~'^'^. We used the peak flux density 
for point sources and the integrated flux density for ex- 
tended sources, and we defined a source as extended if it 
had a ratio of integrated to peak flux density Smt/ Speak 
> 1.1. To account for the varying sensitivity of the sur- 
vey, we also include a correction for the effective survey 
area. The effective area was calculated by weighting each 
source in the bin by the inverse of the area over which 
the source could h ave been detected (|Katgert et al.|[T973l : 
iWhite et al.lll997t ). which we found by summing the area 
in the coverage map with rms < Speak / 5. 

The Stripe 82 source counts up to ~1 mJy are shown 
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Fig. 18. — Stripe 82 source counts below 1 m Jy compared to th e 
source counts from the VLA-COSMOS survey l IBondi et al.|[200l '). 
The counts are shown with and without the correction for varying 
survey sensitivity. 

in Figure [18] compared to the source cou nts from the 
VLA-COSMOS survey (IBondi et all 120081 ). The Stripe 
82 counts are shown with and without the correction for 
varying surv ey sensitivity . The VLA-COSMOS source 
counts from IBondi et al.l ()2008[ ) also include this cor- 
rection, as well as empirical corrections for bandwidth 
smearing, noise bias, fitting errors, and resolution bias. 
The difference between the VLA-COSMOS counts and 
the Stripe 82 counts therefore indicates the remaining 
incompleteness in the Stripe 82 catalog. 

The Stripe 82 survey is clearly incomplete for sources 
near the detection limit. The source counts are a fac- 
tor of 3.5 lower than the VLA-COSMOS source counts 
for the lowest flux density bin ((S) = 0.34 mJy). The 
discrepancy decreases with increasing flux density, with 
the Stripe 82 source counts coming into agreement with 
VLA-COSMOS around S = 0.6 mJy. The drop-off in 
source counts for low flux density bins is most likely due 
primarily to the resolution bias, and catalog users should 
be aware that the catalog is incomplete for faint and ex- 
tended sources. 

9. OPTICAL IDENTIFICATION 

Optical identifications for the Stripe 82 radio sources 
were determined by cross-matching the radio catalog 
with the deep SDSS imaging of Stripe 82 introduced in 
Section HI As this coverage does not extend quite as far 
in Declination from 5 = as our radio coverage, for the 
purposes of this section we trim our catalog so as not to 
falsely underestimate the optical matching rate. 

To determine a suitable matching radius between the 
two catalogs, we first examined the matching complete- 
ness and efficiency as a function of matching radius. For 
this exercise, we used a liberal matching radius of 15" and 
kept the closest SDSS match to each radio source within 
this radius. One-to-one matching was required, with 
SDSS sources flagged once they had been matched once. 
To model the underlying background of false matches, 
we shifted the radio source positions by 1' in four differ- 
ent directions and repeated the matching, averaging the 
results. 

Figure [T9l shows the matches as a function of matching 
radius in bins of 0.3". The total matches are shown by 
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Fig. 19. — Matching our radio catalog to the deep SDSS imaging 
of Stripe 82: Number of matches versus matching radius. The black 
line shows the total number of matches, and the red line shows the 
background matches calculated by shifting the radio catalog by 
1' in four directions and repeating the matching, averaging the 
results. 

the black line, and the background matches are shown by 
the red line. The total matches peak at N = 4816 for r < 
0.3", but the y-axis cuts off at N = 1000 in order to show 
the background more clearly. The difference between 
the total and background matches is then an estimate 
of the number of matches that are actual physical asso- 
ciations. Both the total and background matches show a 
turnover and eventual decrease in matches with increas- 
ing matching radius, a result of the one-to-one matching 
constraint. Modeling of the chance background rate in 
the way described above is not completely accurate, as 
real matches cannot also be false matches. Simply shift- 
ing one of th e catalogs theref o re ove r-predicts the false 
match rate. iMcMahon et al.l ()2002f ) discuss this effect 
in more detail and derive an expression for the effec- 
tive background source density. We have applied a first- 
order approximation of this correction to the background 
matches shown in Figure [T9l where the true background 
is then given by: 

iVw ^ N^.s f ^oni^^'^^^^ - ESbackground \ 

where Nobj are the number of objects in the radio cata- 
log and the summations go up to some radius, r, where 
the background is not yet dominating the total matches. 
We use a radius of r = 1" to calculate the correction, 
and the corrected background and total matches show 
good agreement at large values of the matching radius, 
validating the correction. 

The total matches show a noticeable surplus over the 
background matches at separations beyond that where 
the clear peak in total matches has dropped off. This ex- 
cess, which persists out to separations o f 5" and beyond, 
is not inconsistent with that found by IMcMahon et al.l 
(|2002D for matching between the FIRST radio catalog 
and the SERC Automated Plate Measuring Machine 
(APM) scans of the Palomar Observatory Sky Survey 
(POSS-I) plates, and it likely has two origins. The first 
is that, as radio galaxies often reside in galaxy clusters, 
by allowing matching separations out to 15", we could 
be matching radio sources to unassociated but nearby 
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Fig. 20. — The completeness (solid line) and reliability (dashed 
line) of matching the Stripe 82 radio catalog to the deep SDSS 
imaging as a function of matching radius. 

cluster members. The true background is therefore un- 
derestimated by simply shifting one of the catalogs and 
re-matching, as the actual source density in these in- 
stances is much higher. The second possibility is that 
components of multi-component radio sources are being 
matched to an optical counterpart which does not align 
with any of the radio components. In this case the as- 
sociation is real in the sense that the optical and radio 
sources are related, but without redshift information or 
a smarter algorithm to definitively tie together the mul- 
tiple radio components based on morphology, the match 
rate of the radio sources will be artificially low. 

Choosing a matching radius always comes down to a 
tradeoff between completeness and reliability, with the 
decision ultimately informed by what sort of science one 
hopes to accomplish. The completeness and reliability 
as a function of matching radius are shown in Figure 
[201 Completeness (solid line) gives the percentage of 
the physical matches recovered at that matching radius, 
while reliability (dashed line) gives the percentage of to- 
tal matches that are likely to be physical. The surplus 
in total matches out to fairly large (^^5") radii is clearly 
evident. The completeness is approaching 98% at that 
radius, and the reliability is 85%. On the other hand, the 
matching is extremely reliable for small matching radii. 
For separations less than 1", 99% of the matches are 
likely to be physical, while 97% are physical for separa- 
tions less than 2". The values of completeness for 1" and 
2" are 82% and 88%, respectively. 

Since we wish to provide a reliable estimate of the op- 
tical identifications for the sources in our catalog, we 
will opt for a small matching radius in spite of the low 
completeness. For this matching, we will therefore use 
a matching radius of 1". Using this radius and an area 
common to both catalogs (see above), we find a total op- 
tical match rate of 44.4%. For comparison, the Australia 
Telescope Hubble D eep Field-South (ATHDF-S) survey 
(jHuvnh et al.l l2008| ) reported an optical match rate of 
67% when they compared their radio catalog to observa- 
tions of the HDF-S by the Cerro Tololo Inter- American 
Observatory (CTIO). This imaging data is similar in 
depth to the co-added SDSS imaging (I = 23.5) used 
here, but the radio catalog has a detection limit of rN^50 
/xJy. They also found that the I-band magnitude distri- 
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bution peaked between I = 19 and I = 20 and decreased 
for higher magnitudes. They argued that this was a prop- 
erty of the radio sources themselves and not due to the 
incompleteness in the optical imaging, since the fall-off 
occurred well above the limiting magnitude of 23.5 in I. 
This implies that, while there is surely a significant popu- 
lation of radio sources with optical counterparts too faint 
to be detected here, the majority are not lurking just be- 
low the limiting magnitude. 

One situation where one would not necessarily expect 
there to be an optical counterpart for each radio compo- 
nent is in multiple-component systems such as double- 
lobed radio galaxies. To find out if the unmatched popu- 
lation are more likely to be members of such systems 
than sources with optical counterparts, we perform a 
search around each unmatched source for nearby neigh- 
bors within the Stripe 82 radio catalog. We find that 
~70% of the unmatched sources have another Stripe 82 
radio source within 30". We then performed the same 
search around all of the Stripe 82 radio sources, coming 
up with 80%. We therefore cannot conclude that the un- 
matched sources are unmatched because they are more 
likely to be part of multiple-component systems. They 
are actually slightly less likely to have a nearby neighbor, 
although the difference is small enough that it may not 
be significant. 

The optical counterparts can be further classified into 
either galaxies or stellar sources based on SDSS photom- 
etry. "Stellar sources" refer to sources that are not ex- 
tended. While a small number of stars emit in the ra- 
dio, the majority of these radio sources are more likely 
quasars. The results break down into: 40.6% galaxies, 
3.8% stellar, and 55.6% unmatched. 

Figure [2T] shows the optical identification of all sources 
in the radio catalog broken up based on radio source 
strength. Shown in this Figure are two separate quanti- 
ties: the overall SDSS match fraction, and the fraction 
of those objects that correspond to galaxies. The ra- 
dio sources have been binned by peak fiux density, with 
the data points marking the beginning of each bin. The 
overall SDSS match fraction, designated by the diamond 
points, shows a slight decrease with decreasing source 
strength, from ~60% for sources over 100 mJy to just 
over 30% for sources between 0.3 and 0.4 mJy. The roll- 
off at the faint end could be a real effect, but it is most 
likely due to an increasing population of spurious sources 
at fainter fiux densities. However, the fraction of matches 
increase from about 4 mJy to 0.4 mJy, showing that the 
Stripe 82 catalog is very reliable down to 0.4 mJy. 

Meanwhile, the galaxy fraction, designated by the 
crosses, shows a smooth increase with decreasing source 
strength, from '^65% for galaxies at the bright end 
to '--^95% at the faint end. This is to be expected, 
as it is well-known that the source population sam- 
pled by bright radio surveys is primarily composed of 
quasars and powerful radio galaxies, while fainter sur- 
veys sampling the millijansky level and below are in- 
creasingly dominated by a mixture of low-luminosity 
radio-quiet A GN and star- forming ga l axies (for e.g., 
Afonso et all r2 006: Sm olcic et all 120081 : iSevmour etahl 
20081 llbar et al.. ,2009 '). Nevertheless, we do expect to 
detect previously-unknown radio quasars with this sur- 
vey. To see how many new quasars the survey detects, 
we once again turn to a comparison with the SDSS DR5 
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Fig. 21. — Fraction of Stripe 82 radio sources versus 1.4 GHz 
peak flux density in mJy. The diamonds give the fraction of radio 
sources that have an SDSS match within 1". The crosses give the 
fraction of the matched radio sources that have a galaxy as their 
optical counterpart based on SDSS photometry. The remaining 
sources match to stellar counterparts, the majority of which are 
quasars. 

QSO catalog. 

Recall that the FIRST catalog matches to 229 out 
of 3885 quasars from the SDSS DR7 Quasar Catalog 
(jSchneider et al.ll2010[ ) within a matching radius of 5". 
Of those, the Stripe 82 catalog recovers 223. In addi- 
tion to those 223, the new catalog also has radio sources 
matching 76 quasars not previously detected by FIRST. 
The total fraction of spectroscopic quasars that are ra- 
dio sources to the depth probed here is therefore 7.7% 
(± 0.4%). Conversely, the fraction of radio sources that 
match to spectroscopic quasars is 1.7% (± 0.2%). To de- 
termine the reliability of these matches, we applied the 
same analysis as above to characterize a match between 
our catalog and the SDSS DR7 Quasar Catalog. While 
the statistics were not adequate to produce detailed plots 
of reliability and completeness as a function of matching 
radius, simulating the background counts tells us to ex- 
pect only <1 chance coincidence within the 5" matching 
radius. Therefore, at least 75 of these matches are very 
likely genuine associations. While the majority of the 
newly-detected quasars have sub-mJy peak fiux densi- 
ties, five of them have flux densities above 1 mJy, and 
therefore should have been seen by the FIRST survey. 
In three of the cases, however, the source has a double- 
lobed morphology, and FIRST'S resolution was too poor 
to resolve the core from one of the radio lobes. The po- 
sition of the FIRST source was therefore shifted further 
than 5" from the optical center. In the other two cases, 
the FIRST cutout server shows that the sources were 
clearly not bright enough to appear in the FIRST cat- 
alog, implying that the peak fiux density of the quasar 
has brightened in the high resolution data due to variabil- 
ity. All 76 quasars are therefore newly-discovered radio 
quasars, and a list of the sources is presented in Table [2j 

As a further comparison, we also matched our radio 
source catalog to a larger catalog of photo metrically- 
select ed quasars based on the SDSS DR6 ( Ric hards et al.l 
l2009f ). There are ~13,500 photometrically-selected 
quasars in the area covered by our data, almost 3.5 times 
the number of spectroscopically-confirmed quasars. De- 
spite this fact, the fraction of radio sources that match 
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TABLE 2 

Newly Detected Radio Quasars 



RA° Dec Spk' Si„i'= Maj-* Min'' zf mr<i 
(mJy) (mJy) (") (") (mag) 



''NOTES. - "Radio source position {J2000). ^1.4 GHz peak 
flux density (mJy). '^1.4 GHz integrated flux density (mJy). 
"^Deconvolved major axis. "D econvolved minor axi s. -^Redshift 
from the DR7 quasar catalog ISchneider et al]|20ig) . ^ Apparent 
magnitude in the r-band from the same catalog. 
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Fig. 22. — Fraction of quasars that are radio sources as a func- 
tion of SDSS r-band magnitude. The solid line represents the 
spectroscopically-confirmed quasars, while the dashed line shows 
the photometrically-selected quasars. Error bars were calculated 
using the formal binomial distribution: = [/(I — /)/Af]'''^. 

to quasars only increases to 2.0% (± 0.1%), a factor of 
1.2. Figure [22] shows the fraction of quasars that are ra- 
dio sources as a function of SDSS r-band magnitude for 
both the spectroscopic and photometric samples. The 
two samples are consistent above r ~ 20, lending cre- 
dence to the reliability of the photometrically-selected 
sample in that magnitude range. Below r = 20, we see 
a statistically-significant deficit in the radio fraction of 
the photometric sample with respect to the spectroscopic 
sample. As the peak of the photometrically-selected 
quasar magnitude distribution lies in this range, this 
seemingly small difference explains the large discrepancy 
in the total percentages referenced above. 

We believe that this difference can be explained, at 
least in part, by the efficiency of the SDSS target se- 
lection algorithms. Targets may be selected for a num- 
ber of reasons, including stellar objects within 2" of a 
FIRST source, ROSAT sources with quasar colors or that 
are also radio sources, and, perhaps most notably, in an 
open category of serendipitously selected targets includ- 
ing objects coincident with FIRST sources but fainter 
than the equivalent in quasar target selection (and not 
restricted to point sources). As a result, the photometric 
catalog would mainly fill in those objects without ra- 
dio detections, which are presumably much more numer- 
ous. This would explain why the total percentage of 
photometrically-selected quasars that are radio sources 
(based on our catalog) is only 2.7% (± 0.01%). Further 
investigation into the radio and optical properties of the 
entire quasar sample will be undertaken in an upcoming 
paper. 

10. PUBLIC ACCESSIBILITY 

All images are available in searchable format via a link 
on the FIRST websit^B- The search parameters required 
are the same as those required by FIRST and are de- 
scribed online. The resultant cutout will be shown along- 
side the corresponding cutout from the FIRST survey. 
The Stripe 82 catalog is also available online in a search- 
able format0. In addition, an option has been added to 

* http:/ /third. ucllnl.org/cgi-bin/stripe82cutout 
^ http:/ /sundog. stsci.edu/cgi-bin/scarchstripe82 
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the FIRST catalog search page to view the correspond- 
ing Stripe 82 image in the two non-contiguous areas of 
coverage. The full Stripe 82 catalog will be available for 
download shortly from the same website. 
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